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SUMMARY 

Geomagnetic field lines are traced into space for more than 
500 geophysical stations, using the 48-term Jensen and Cain field 
expansion for 1960. Total field intensity, magnetic dip, declination, 
and L-values are computed at each station, and at points 300 and 
400 Km vertically above. Coordinates and B values of the conjugate 
intersects and of the equatorial (minimum B) points are determined 
for each field line, as well as the total arc length of the line between 
origin and conjugate intersect. 

A general discussion of the resulting tables is given. 

Finally, a series of maps showing constant B and L contours 
at various levels between surface and 3,000 Km are reproduced. 
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COSPAR's Panel on Conjugate Point Experiments was established in 1962 
as a subgroup of Working Group 2, on the basis of a recommendation of the 
IQSY Working Group on Cosmic Rays and Geomagnetic ally Trapped Particles. 
Members of the Panel are: R. A. Helliwell, W. N. Hess, J. D. Kalinin, J. Paton, 
R. Schlich, R. L. O. Storey, and E. Vestine. 

So far, the Panel's activities were mainly limited to keeping COSPAR in- 
formed about conjugate point experiments under course, to propose pertinent 
recommendations on cooperative programs and to discuss possibilities for future 
experiments. 

At last, COSPAR's Panel on Conjugate Point Experiments has started to 
produce information for more general use. This was possible thanks to the most 
efficient cooperation of NASA's Goddard Space Flight Center, which is herewith 
duly acknowledged. 

We believe the following tables may find wide use and may constitute the 
beginning of our Panel's activities as a service to the scientific community. 
Moreover, it is our sincere hope that these tables might serve to increment 
cooperation or to start new scientific ties between different countries — in this 
case "magnetically conjugate countries" a service which, after all, is the 
ultimate goal of all COSPAR Working Groups and Panels. 


Juan G. Roederer 
Chairman 

COSPAR Panel on Conjugate Point Experiments 
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CONJUGATE INTERSECTS TO 
SELECTED GEOPHYSICAL STATIONS 
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INTRODUCTION 

In the past few years, the study of propagation of charged particles, elec- 
tromagnetic radiation and magnetohydrodynamic waves along geomagnetic field 
lines has become a powerful tool in the search for a better knowledge of the 
structure of the Earth's magnetosphere. 

The determination of the exact shape of a line of force, its change with local 
time and during geophysical perturbations, and the answer to the question of 
return to Earth of high latitude field lines, are of crucial importance to the prob- 
lem of the interaction of the solar wind with the Earth's magnetic field. The 
process of propagation of magnetohydrodynamic and very low frequency electro- 
magnetic waves along field lines gives information on the electron and ion den- 
sities encountered; the study of precipitation of energetic particles down the field 
lines into the atmosphere hopefully will lead to a better understanding of the 
injection and acceleration mechanisms in the magnetosphere. 

Many of these investigations are based on ground measurements. Hundreds 
of IQSY and ex-IGY stations gather information which is being used or which 
could be used for these studies. These measurements are particularly valuable, 
if carried out simultaneously at both "ends" of the same field line. It is there- 
fore very useful to know with maximum possible precision, the geometry of the 
field line which "goes through" each station, or which intersects the ionosphere 
at several given vertical altitudes, for a given field expansion and during normal, 
quiet geomagnetic conditions. 

There are many useful parameters related to this "field line geometry." 

For this preliminary edition, we have selected the following: 

1. Geographic coordinates of the station (or the vertical point at high 
altitude), henceforth called the origin of the field line; 

2. Total field intensity at origin; 

3. The magnetic dip at origin; 

4. The magnetic declination at origin; 
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5. L -value corresponding to the origin (and within about 1%, to the whole 
line); 

6. Geographic coordinates of the conjugate intersect, i.e. the intersection 
of the field line with the Earth's surface (or a fixed level at the same 
altitude as the origin), at the other side; 

7. Total field intensity at the conjugate intersect; 


8. Geographic coordinates of the equatorial point of the field line (defined 
as the point where the field intensity attains its minimum value on that 
line); 

9. Total field intensity at the equatorial point of the field line; 


10. Length of the field line between origin and conjugate intersect. 


These quantities, which are printed out in the above order in Tables 2,3, 
and 4, by no means represent all useful field magnitudes. However, for reasons 
of space and time, other important quantities are not reproduced here. 

Stations included in our tables were selected from IQSY lists, as well as 
extracted more or less at random from old IGY tables. This selection by no 
means pretends to be complete; we apologize for omissions and possible errors 
in name -spelling, coordinates, etc. 

A careful analysis of the geographic location of conjugate areas was made. 
About 380 stations have their conjugates over sea. Only 140 stations have con- 
jugate intersects over land, 16 of them over Antarctica. And only 22 pairs of 
stations are nearly self -conjugate, some of them deliberately established as such, 
others just by chance. We present a list of these self-conjugate pairs, and give 
a survey of the geographic location of all other conjugate intersects which fall 
over land, indicating country, and whenever possible, the nearest town to the 
intersect. 

Finally, for the sake of completeness, we present a series of maps with 
constant B and L contours, for altitudes between sea level and 3,000 Km. 


FIELD LINE PARAMETERS 


1. Selection of Stations 


All Stations were selected from lists published in IQSY Notes No. 4 (1963) 
and No. 11 (1965), as well as extracted more or less at random from the IGY 
Annals. Some stations of this latter group may be at present out of operation. 

Whenever a group of stations appeared lumped together in the same geo- 
graphical area, only one or two of them were picked out (generally the station 
bearing the name of the near by greater city — example: Of the pair of stations 
Godley Head and Christchurch, only the latter was kept). 

For practical reasons, stations are given in our lists in alphabetical order, 
according to the first two letters. 


2. Coordinates 

All positions are given in geographic coordinates . Latitudes and longitudes 
are in degrees and decimal fractions; negative latitudes are South, negative 
longitudes are West, Altitudes and arc lengths are given in kilometers. 

Information on altitude over sea level was not available for many stations. 

In case of doubt, sea level was adopted. A test has shown that a change of a few 
thousand meters in station altitude does not alter appreciably the values listed 
in our tables, except for stations at very low geomagnetic latitudes. 

All calculations were performed for field lines passing through the stations, 
and through points 300 and 400 Km vertically above each. These two levels, 
which represent typical heights of the F2 layer in its diurnal variation, were be- 
lieved to be of interest to ionospheric researchers. They may serve to obtain 
interpolated values for field line parameters of lines of force passing through 
other heights than the ones selected here. 


3. Field Parameters 


The geomagnetic field used in these calculations is the 48 term expansion 
by Jensen and Cain 1960.0. 1 A more recent field expansion by Cain et al., 2 
which takes into account the oblateness of the Earth and updates the field to 1965, 
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was applied to several stations for comparison.* Differences were very small; 
conjugate points fell within 20-30 Km of the positions given in our tables. 

For the actual computation of the field, Mcllwain's subroutine MAGNET 
was used.** Values of the field are given in gauss/ The error is expected to 
be of the order of 2 x 10 -3 gauss. 

The conjugate intersect is defined as the intersection of the field line emerg- 
ing from the origin (surface, 300 Km or 400 Km above s.l.), with a constant level 
of the same height as the origin. (See Figure 1). 

The "minimum-B equator" is defined as the point of minimum total field in- 
tensity on a given line of force. The locus of all these points is, of course, not 
coincident with the magnetic dip equator. 

The altitude of the minimum-B equator is very nearly equal to the maximum 
altitude of the field line except for low L lines (L < 1.2). Their latitudes, how- 
ever, may differ appreciably. This was tested in separate runs for all stations. 

The longitude of the minimum-B equator is the parameter "par excellence" 
to label a given field line on a L shell, in a nearly invariant way. 4 

For very low geomagnetic latitudes, it sometimes happens that there is no 
minimum-B value between origin and conjugate intersect. This means that the 
minimum-B equator of that particular field line lies below the origin, or below 
the conjugate intersect (Figure 2). This may even occur in the case of surface 
to surface tracing; it is an indication that in the region considered, there is a 
difference between the "magnetic" and the "minimum-B," or "invariant," equator. 

Finally, the L-parameter, 3 given in Earth radii, is computed through 
Mcllwain's INVAR subroutine.* 


*We are very much indebted to Dr. J. Cain for carrying out these calculations. 

**MAGNET is port of INVAR (see footnote below). 
tAs we use “B”, we should call the unit Oersted. 

t INVAR is a computer code written by C. E. Mcl twain, which has been widely distributed. 
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If L is high, say, above 8 earth radii, the field liqe expansion used in MAG- 
NET is no longer valid, because external sources, as determined by the inter- 
action of the solar wind with the earth's magnetosphere, can then no longer be 
neglected. This is the reason why field lines were traced only when L £ 12. In 
general, above L = 8, this parameter merely has to be taken as a quantity re- 
lated to the "invariant" latitude defined by cos 2 A = l/L. 

It is well known 3,5 that L is constant for a given magnetic shell only within 
about 1%. Azimuthal asymmetry of the real geomagnetic field actually splits 
particle shells according to their equatorial pitch angles (or mirror point field 
intensities) at a given field line. The INVAR-computed L-values fluctuate on a 
given "L-shell," and even along a given field line. A good measure of this fluc- 
tuation is the quantity (L-Lo)/Lo where Lo = (Mo/Bo) 1/3 (M = 0.3117, Bo = min- 
imum field intensity on the given field line). A test made with our results has 
shown that the longitude dependence of (L-Lo)/Lois in excellent agreement with 
Mcllwain's results (Figure 3 in 3 ). A significant difference, however, arises in 
the region of longitudes between 5 and 35 degrees East, for L values in the in- 
terval 1.6 to 2.6. There we find computed L values 2 -4% higher than the "azi- 
muthally symmetric" value Lo. 

For many purposes, it may be sufficient approximation to use Lo for the L 
value of a station, or of a point in space. The elimination of the subroutine 
INVAR from the program would then save considerable computer time. 

In order to make our tables useful to cosmic ray physicists, we may point 
out the results of Smart and Shea, 6 who found an excellent correlation between 
L and the effective geomagnetic cosmic ray cutoff rigidity IJ. . We have reana- 
lyzed this correlation using the computed surface L values given in Table 2, and 
concluded that the effective cutoff rigidity can be found for each station by the 
expression 

P c = 15.60 L" 2 (GV) (1) 

valid for L £ 1.2. For equatorial regions this correlation breaks down. 


METHOD 

1. Coordinate System Used 

Figure 3 illustrates the system used. The field components along the unit 
vectors u v , u N and u £ are: 


B v (vertical downw ards) 
B n (horizontal North ) 

B £ (horizontal East ) 
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Figure 3 


The equations defining an infinitesimal portion of a line of force are then: 

SS _ _Sr_ _ r S A. _ rcos^Sy 

B By B n B e 


SS is the element of arc, with components — Sr, rSA. and rcosXScp along the 
unit vectors. The coordinates of a generic point of a field line originating at 
r s , \ and cp s are then given by 


r 


= r 

n s 



(3a) 


k = k. 



— SS. 
B. 


9n= 9 



J 

COS k . 

1 



SS. 


(3b) 


(3c) 


where the sign factor cr is ±1 according to whether one proceeds tracing in the 
opposite/same direction of the field vector. Notice that if the origin is in the 
magnetic Northern hemisphere (By > 0 ), we have to set o' > 0 in order to follow 
the line to the conjugate intersect* (direction opposite to B ). 


*ln our fieid computations using subroutine MAGNET, we obtain the West component of the field 
B w =-B £ . For this reason we had to reverse the sign in (3c). (See listing. Table 1). 
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The total arc length between origin and a generic point is 


% 



(4) 


COMPUTER PROGRAM 
1. Procedure 


The flow diagram is shown in Figure 4. 

The altitude level of the origin is determined by an input parameter which 
permits one of three internally programmed choices (station altitude, 300 and 
400 Km). Control parameters are initialized and the station is entered. 

The L value is computed with Mcllwain's INVAR subroutine. The L value 
is then tested; if L > 12, no tracing is attempted. Before the program returns 
to read in the next station, it computes the total field, the magnetic dip and the 
declination at the origin. 

If L < 12, the code proceeds to trace the field line. The "Adams Four Point 
Integration Formula" 7 for numerical solutions of differential equations is used 
in this procedure repeatedly at equidistant intervals. The accuracy of this fit 
depends on the integration step size (element of arc length). 

To increase efficiency and to save computer time, the constant step size D 
for a given field line was fixed as a linear function of L: 

D = 50 L - 30 (Km) (5) 


Exhaustive testing with several step sizes for fixed L values has shown that 
relation (5) is satisfactory with respect to both accuracy and machine time. 

Initial conditions for the integration of the dependent variables H (altitude), 
DLAT and DLON are given, and the total field and its components are computed. 
The dip and the declination are calculated only once, on the first passage through 
the integration loop. A suitable control causes this part of the program to be 
bypassed thereafter. 
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2. Usage and Comments 


The computer code is a FORTRAN II monitor controlled program written 
for the IBM 7090-7094, Standard or Moonlight system (Table 1). 

Execution time varies from fractions of a second for low latitude stations 
(small L-values) to several seconds for stations at high latitudes (large L-values). 
For the 511 stations considered in this work, the average time was approximately 
2 seconds per station. 

The choice of altitude levels presently contained in the code is explained on 
page 3. Should other considerations prevail or different use be made of the 
routine, new values can easily be substituted in statements number 58 and 59. 

An effort has been made to keep the code flexible for diverse usage. Addi- 
tional quantities connected with a field line, such as: 

a. maximum altitude attained above the earth's surface, 

b. geographic coordinates of the intersect with the geographic equator, 

c. conjugate mirror-points, 

d. curvature of the equator, 

e. field components at origin and conjugate intersect 
may become readily available through minor modifications. 

The field computing part of the program has been kept separate in the form 
of a subroutine. This is practical and advisable, as experience has shown. Since 
there exist several models describing the geomagnetic field in terms of a varying 
number of expansion coefficients and since periodic recomputations of these 
coefficients are necessary due to a gradually better knowledge of the field and 
its secular changes, this arrangement is most convenient for interchanging or 
updating the various models. 

Two cases may arise which receive special treatment in the routine: 

a. The origin may have an L-value greater than 12, 

b. No minimum-B equator may lie between the origin and the conjugate 
intersect. 
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In both, the code will indicate their occurrence by printing relevant comments 
in place of the regular output. 

In the tracing process, the location of the minimum-B equator is approxi- 
mated on the field line with an error proportional to plus or minus one integration 
step. 

There are no limitations with respect to the number of stations or the height 
of the origin. 

The restriction L 1 12 imposed on the shell parameter is due to physical 
(see page 5 ) and practical considerations only (to save computer time) and is 
entirely independent of the code itself. Lines with L-values greater than 100 
have been traced repeatedly without any difficulty. 

In the loop, the following sequence of operations is performed: (a) the total 
field intensity is tested at every step, and the coordinates and field intensity of 
the minimum-B point are stored; (b) the altitude, latitude and longitude increment 
parameters are evaluated; (c) the n-1 th dependent and independent variables are 
stored; (d) the nth numerical integration is performed. 

The integration is terminated whenever the altitude H n is less or equal to 
the altitude of the origin H s . The total arc length is then tested; if S < 1000 Km, 
the step size is decreased by a factor of 4, and the tracing is repeated for that 
given field line. This is done in order to obtain a better fit and to reduce the 
error in the location of the conjugate intersect, for low latitude stations. 

If S > 1000, or after a second tracing, the coordinates of the conjugate 
intersect are obtained by interpolation, and the field intensity is computed at 
that point. 

Finally, a test is made to determine whether the minimum-B equator lies 
between the origin and the conjugate intersect. The results are printed, and the 
next station is read in. 


RESULTS 

Table 2 gives the results of our computations for surface-to-surface tracing 
for 511 geophysical stations. In column 1, the names of the stations are listed. 
Columns 2 and 3 give the geographical coordinates of the origin; column 4, the 
total field in gauss; columns 5 and 6, magnetic dip and declination in degrees 
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respectively; and column 7, the L value. Columns 8 and 9 give the coordinates 
of the conjugate intersect, column 10, the total field at the intersect. Columns 
11 and 12 show the coordinates of the equatorial (minimum B) point of the corres- 
ponding field line, column 13, the field intensity there, and column 14, the altitude 
of this point over sea level (almost coincident with the maximum altitude over 
sea level of that line, except for low L values). Column 15, finally, gives the 
total arc length of the field line, between origin and conjugate intersect. All dis- 
tances are in kilometers. 

Tables 3 and 4 reproduce the values corresponding to field lines originating 
300 and 400 Km above each station and intersecting the same altitude as the 
origin, at the conjugate point. 

A detailed inspection of the conjugate areas has revealed several features. 
First of all, as expected, the great majority of the stations have conjugate areas 
over the sea. There are, nevertheless, 22 pairs of stations which are very nearly 
conjugate to each other. They are listed in Table 5. For additional information, 
the mean L value is given, as well as the percentage difference in L. With a 
few exceptions, this difference is not much greater than the intrinsic fluctuation 
of L. Another important parameter necessary to judge the "degree of conjugacy" 
of two stations, is given by the difference of longitudes of the equatorial points 
of the field lines passing through both stations. This is indicated in the last 
column of Table 5. 

Inspecting Table 5 one recognizes well-known conjugate pairs, in addition 
to others. Notice in particular, the European-African group of conjugate pairs. 

Table 6 lists stations with conjugate areas in the Antarctica. Those already 
listed in Table 5 are not shown here. 

Finally, in Table 7 we give some information about the conjugate area for 
those stations whose conjugate intersects fall over land, but where there is no 
other station near by, belonging to our list. In this table, we list the name of 
the station in alphabetical order, the L value, the name of the city or of the region 
in the conjugate area, and the country in question. 


B-L CONTOURS 

The last collection of figures represent constant B and L contours for differ- 
ent altitudes. They were obtained in the following way: For a given altitude 
level, a latitude -longitude grid is generated with constant intervals of five de- 
grees for the first and of ten degrees for the latter. At those grid points, the 
magnetic parameters B and L are computed and stored. 
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The desired constant B or L contours are then obtained by a linear inter- 
polation, which is performed separately for both. The resulting positions are 
then plotted by computer and the curves are drawn manually. 
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Table 1 

Arguments and Parameters 




I. Input: 


J 


BLAT 

BLON 

BH 

NAME 


II. Internal : 

RAD 

N 

LP 

ITEST 

D 

S 

H 

DLAT 

DLON 

SGN 

DECL 

HM 

B 

BN 

BW 

BV 

FAC 

DH 

DDLAT 

DDLON 

CBN 

CBW 

CBV 


fixed point parameter controlling starting altitude; for 

J = 1, 2, 3 altitude is fixed in the code to station height, 

300, 400 Km, respectively. 

geographic latitude of station. 

geographic longitude of station. 

height of station above sea level. 

designation of station. 


conversion factor, radians to degrees. 

fixed point parameter counting lines per output page. 

fixed point parameter numbering tracings per line. 

fixed point parameter controlling dip angle and declination 

bypass 

constant incremental change of independent variable 
(stepsize of integration along field line), 
independent variable: total arc-length, 
dependent variable: height above sea level, 
dependent variable: geographic latitude, 
dependent variable: geographic longitude, 
sign of vertical field component, 
magnetic declination at origin. 

altitude above sea level of point of minimum B on field 
line. 

field strength at generic point. 

North component of field vector at generic point. 

West component of field vector at generic point 
(See footnote page 6). 

vertical component of field vector at generic point. 

factor for increment parameters DDLAT and DDLON. 

altitude increment parameter. 

latitude increment parameter. 

longitude increment parameter. 

same as BN at conjugate intersect. 

same as BW at conjugate intersect. 

same as BV at conjugate intersect. 
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Table 1 (Continued) 


♦ 


SP 

DLATP 

DLONP 

HP 

C2 


^u-xjui reference arclength for conjugate intersect 
approximation 

(n-l)th reference latitude for conjugate intersect 
approximation 

(n-l)th reference longitude for conjugate intersect 
approximation 

(n-l)th reference altitude for conjugate intersect 
approximation 

ratio of B(min) minus B(orig) to B(orig). 


HI. Output: 


BB 

DIP 

IDECL 

BEL 

CL AT 

CLON 

CB 

DLATM 

DLONM 

IHM 

BX 

DS 


field strength at origin, 
dip-angle at origin, 
see DECL 

magnetic shell parameter L of field line through origin, 
geographic latitude of conjugate intersect, 
geographic longitude of conjugate intersect, 
field strength at conjugate intersect, 
geographic latitude of minimum B point on field line, 
geographic longitude of minimum B point on field line, 
see HM 

minimum field strength on field line. 

total arclength between origin and conjugate intersect. 


NOTE : All arguments are in floating point form, except where otherwise 
indicated. 
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Table 1 

FORTRAN LISTING 


* CARDS COLUMN 

* LI 01 8 

* LABEL 
CL0F1A 

MODIFIED LuF TRACER » VERSION NO. 3 FOR ALL STARTING ALTITUDES 

FIELD COMPUTATION BY MCILwAINS MAG ME 1 

PARAMETER J HAS FOLLOWING VALUES AMD MEANING = 

J=1 LOF TRACING FROM SURFACE TO SURFACE 

j=? lof tracing 3 uo km abovf station to 300 km conj. intersect 

jr 3 LOF tracing AUO km above station to aoo km con j. i ntersect 

DIMENSION NAMFC*) 

READ I NPU I Tape 2*1uua*J 
1UUA FORMAT (13) 

RAD=57. 2957795 

K = 1 
N = 0 

1 READ INPUT TAPE 2*lo05*BLAl* BLON »BH » NAME 
1UU5 FORMAT ( 2F 1 1 . 2 * F 8 . 2 * 3 A A ’ 

I TEST = 1 
LP = 1 

IE ( J - 2)57*58*59 

58 bH = 30u . 

GO TO 57 

59 BH = AGO. 

57 CALL INVAR(BLAT*BLON*BH *0.01 * BE> > BEL ) 

IF (BEL - 12.12*2*151 
151 SENSE LIGHT A 
GO TO 6 

2 D = 50 • *BEL - 30. 

6 DLAT = BLAT 
DLON = BLON 
H = BH 

91 S = 0 .0 

SENSE LIGHT 1 
BX = 5. 

92 CALL F I ELD ( DLAT *DLON *H» BN *BW *BV *B ) 

GO TO (38*39)* I TEST 

38 DIP = ATANF (BV/SQRTF (BN**2 + BW**2))*RAD 
uGN = S I GNF ( 1. *BV ) 

DECL =(ATANF(-BW/BN))*RAD 
IDLCL = DECL + .5 
37 ITEST = 2 

I F ( SENSE LIGHT A) 7*39 

7 IF(N)9*9*1U 

9 WRITE OUTPUT TAPE 3*73 
WRITE OUTPUT TAPE 3*7A 

lu WRITE OUTPUT TAPE 3 * 77 * NAME ♦ BLAT » BLON *BB » D I P * I DECL *BE L 
77 FORMATl 1X3A6*F8.2 *F9.2*F 7.3*E6. 1 * I 6 *F7.2 *AX55H** NO C 0 M P U T 
1 A T IONS PERFORMED**) 

N = N + 1 
GO TO 1 

39 I F ( B - B X ) 3 » A * A 

3 DLATM=DLAT 
DLONM=DLON 
HM = II 

BX = B 

A DH = ( BV/B ) *SGN 

FAC = SIGNF (RAD/ (B* (11+6371 .2 ) ) *SGN ) 

DDLAT =-bn*fac 

DDLON= BW*F AC/COSF < DLAT /RAD ) 
bp = b 
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** HP = H 

DLATP = DLAT 
DLONP = DLOM 
b=INDVF (b.D) 

H=DPNVF ( H *DH ) 

DLAT =DPNVF( DLAT, DDLAT ) 

DLON=DPNVF ( DLON . DDLOri ) 
i F ( i I - BH ! 22 * 2 2 * 92 

22 I F ( S - 1000. ) 2 6 > 2 6 » 5 
26 GO TO ( 150.5) »LP 

150 D = D/4. 

LP = 2 
GO TO 6 

5 DLAT = DLATP + (DLAT - DLATP )*<BH - HP)/(H - HP) 

DLON = DLONP + ( DLON - DLONP)*(BH - HP)/(H - HP) 

Db = bP + (b - bP ) * ( BH - HP ) / ( H - HP) 

H = BH 
CLAT = DLAT 
CLON = DLON 

CALL FI ELD (CLAT .CLON . H , CBN > CBW , CBV » CB ) 

IHM = HM + .5 
IF (CB - BX >51.51.160 
16b C2 = BX/BB - 1. 

52 I F ( AfabF ( C 2 ) - 1 . E-5 ) 5 1 . 5 1 » 20 
51 bENbE LIGHT 2 
bENbf LIGHT 3 

2b I F ( ABSF (CLON) - 180 . ) 88 . 88 . 87 

87 I F ( CLON >82.88.83 

82 CLON = CLON + 360. 

GO TO 88 

83 CLON = CLON - 360. 

88 IFlbENbF LIGHT 3)46.45 

45 I F ( ABbF ( DLONP) ) - 1 80 . ) A6 .46 . 84 

84 IF(DLONM) 85*46.86 

85 DLONM = DLONP; + 360. 

GO TO 46 

86 DLONM = DLONP) - 360. 

46 IF(N)3o, 30.21 

3w w R IT E OOTPOT TAPE 3»73 

73 FORMAT (131H1 bT AT I ONb ** 0 R I G I 

1 N ** CONJUGATE I NT ERG . ** i'ilN IMUM-B EQUATOR ** AR 
2C //) 

WRITE OUTPUT TAPE 3.74 

74 FORMAT ( 22X4HLAT. . 5X5HLONG. » 4X 1HB . 4X3HD I P » 3X5HDECL . . 3X 1 HL . 5X4HLAT . , 
13X5HLONG. .4X1HB.5X4HLAT . .4X5HLONG. .4X4HALT . » 4X II IB * 6X2HKM/ / / ) 

21 I F ( bENbE LIGHT 2)23*31 

23 WRITE OUTPUT TAPE 3 * 78 » NAME » E-.LAT .PLON »BB »D IP . I DECL .BEL .CLAT .CLON . C 
lb.Db 

78 FORMAT l 1X3A6.F 8.2.F9.2. F7.3»F6.1»I6»F7.2.F7.1»F8.1 . F 7 . 3 » 4X28HNO MI 
INI PiUP'i b BETwELN I NT Ekb. » F8.0 ) 

GO TO 90 

31 WRITE OUTPUT TAPE 3» 75 » NAME » BLAT . BLON .BB.DIP. IDECL.BLL.CLAT * CLON . C 
1 B . DLATP) . DLONM . I HM . BX , Db 

75 FORMAT (1X3A6.F8,2*F9.2»F7.3.F6.1.I6*F7.2*F7.1*F8,1»F7.3»F8.2»F9.2» 
1I7.F7.3.F9.0) 

9u N = N + 1 

I F ( N - 45)32.33.33 
33 N = 0 

WRITE OUTPUT TAPE 3.76.K 

76 FORMAT ( ////6uX5HPAGL 14) 


K = K + 1 
32 GO TO 1 
END 
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Table 2 

Tracings from Surface to Surface 
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Table 2 

Tracings from Surface to Surface 
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Tracings from Surface to Surface 
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Tracings from Surface to Surface 
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Table 2 

Tracings from Surface to Surface 

STATIONS *• 0 R I G I N ** CONJUGATE INTERS, ** MINIMUM-B EQUATOR ** ARC 
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Tracings from Surface to Surface 

STATIONS •• o k l G I N *• CONJUGATE INTERS, •• HINIMUM-B EQUATOR •« ARC 
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PORT STANLEY 


Table 2 

Tracings from Surface to Surface 


U 

an 

< 


o 

cr 


o 

CO Z 

1 <-J 


x 

z 

a: 


o 

z 


< 

o 

-> 

z 

o 

o 


cr ao x cm (\j 4- 4*X(xj 3— i4"xxoo4'cr aoooxo3aoao30'xx3xxxaoooxa0'. 

Q04'CT4-XCra030’-<Crx(NjXXX3 * rf'i — imr-ioOOO>OMroN®(\jH(\j>OOOirco®( 
(M 0 ' 0 iri<oin 0 si-'j , (\iri>oc 0 r-< 0 '-o^ * O'Om^^oocrinooivjaj^iiuASHO'Oi 

4 - 3 (nj(njxox 4 ’ 04 >' 0 (njx oo o x 3 O' x m 3 x ao 03 X 4 ‘XX 303 Xxco< 
x<ni(\j 33 X 3 mh x in 3 o O' 3 h rn ^ — » vO ^ h ^ h h . 


3 3 3 O' 4* 

4 > in x >o O' oo 

* x in x 4- nj 

<r cr x 3 cr 

a x (xi 3 

LU 


i in 4 

in 

X 

m >j- 

<XJ f\l 

0 

X X X l/l 

0 

r- 0 


_* 

<r 

O' 


X 

(X) 

(XI 

i/> 

X 

3 X 

4- 

4 - x 3 

3 X 

(XJ 

X 

3 

3 X 4) 4" 3 3 

• 0 x 

x 

3 

0 — • 

pM rn 

r~4 

0 >4- 0 CC 

(XI 

O (XJ 

X 

0 

<4- 

1 

00 

r- < 

CD 

X 

CC 

0 

x cr 

O 

0 O' X 

O O' 

(XJ 

rg 

X 

4 - 3l 4- 4- 3 4- X 

i 0 0 

0 

0 

0 0 

O rM 

(XJ 

O X 0 UJ 

H 

3 (Xj 


0 

0 

ru 

0 

O 

O 

0 

3 

*— 1 

3-4 3-4 

O 

000 

3 O 

(XJ 

(XI 

3 

O 3 3 O O O 

> • • 

♦ 

• 

• • 

« « 

« 

• • « h 

* 

• * 

cc 

• 

• 

•- 

• 

♦ 

« 

• 

3 

• 

• • 

• 

• t • 

• • 

• 

* 

♦ 

• cc • • ♦ « • 

0 

0 

0 

0 

0 

0 

0 

0 

O 

0 

0 

0 

N 

0 

0 

0 


0 

0 

0 

0 

O 

O 

0 

Z 

0 

0 

0 

O 

000 

O O 

0 

0 

O 

0 00000 










0 

















a 

> nj 3 

CD 

3 

m ru 

0 r- 

>T 

O s4- Z 

X 

x 0 

3 

0 

nj 

3 

X 

oo 

O 

(XJ 

z 

3 

x x 

X 

3 3 x 

X O' 

3 

X 

3 

X LL CT (XI 0 X X 

j >0 — ' 

in 

x 

>r 0 

r- <3 

X 

(\l 00 P* UJ 

>4- 

cr x 



X 

X 

00 

X 


CO 

3 

4" 

3 X 

CO 

0 3 3 

X X 

O' 

X 


X 4> 3 (M O 3 

) 3 3 

X 

x 

O' r- 

0 in 

fXJ 

X (XJ X 3 

O 

x <r 

CC 


X 

X 

(XJ 

0 

4- 

«— 1 

3 

X 

4" 3 

3 

XXX 

X-3 

(NJ 

(X) 

O 

X CC (XJ (XJ 45 X X 

O' X 

>0 

0 

re — » 

0 rxi 

r— 4 

X X X 

(XI 

(XJ 


re 

X 


X 

00 

X 

X 

X 

(XI 

3 

3-4 

3X4- 

(XJ (XJ 



(XI 

X (XJ (XJ 3 4) X 




(XJ r-4 



-* 3 3 



3 

X 




3-4 



3 



(XI 

(XJ 





UJ 3 







3 











3 
















X 



a 








0Q 









Ol 

r x co 

0 

(XJ 

sD O' 

O' O' 

X 

>4-0 0 

X 

n4- X 


X 

X 

ao 

X 

0 

X 

X 


3 

X (NJ 

(XI 

X X 0 

cr 0 

O 

X 

X 

3 X O' 4> X X 

) aj in 

X 

0 

>o r- 

m in 

>0 

3 cr cr x 

NT 

X (XJ 


X 


QO 

X 

X 

CD 

3 

CD 

X 

X 3 

3 

4* X O 

3 O' 

X 

-j 

O' 

O X 3 X (XI (XI 

3 3 0 

0 


® cr 

OOP- 

so 0 x a: 

0 

vT X 

X 

O 

X 

X 

X 

3 

X 

H 

2- 

O 

O' cr 

O 

cr nj x 

X 4- 

X 

X 

3 

4" X O 3 X 3 X 

- 3 rxi 

(XI 

(XI 

r* (XI 

VU fXJ >0 

(Xj (XJ -o 3 

1 

X X 




>4- 

X 

X 


4- 

3 


1 x 

3 

3 (XI 3 

X X 

4- 

4- 

(XI 

3 4 - (XI X X 3 




^4 -1 

1 


1 — 13 : 


1 -4 

z 



1 

1 




2. 


l 


333 

1 1 

I 

i 


Z 3 3 3 




1 1 



z 



a 





1 



Z 




1 1 





O 




























3 3 3 

(XI 

4 

0 — * 

P- sj- 

X 

O 3 X 

X 

3 X 

km 

O 

QO 

(XJ 

oo 

<4- 

X 

>T 

£ 

X 

(\j ao 

3 

O' X X 

O O 

X 

X 

3 

0 — X 4) 3 3 3 

' (\j O 

X 

CXj 

in *-■ 



X 3 <5 


O (XI 


O' 

0 

X 

X 

■4" 

00 0 


cr 

3 0 

3 

3 m 3 

cr cr 

X 

X 

X 

(XJ 0 O' X 3 4” 


• 

• • 


• 

• • • O 


• « 

3 

• 


• 


• 



a 





• • 




• 3 • • • * • 

J X X 

O 

in 

r- in 

-4- O' 

r- 

P- X fXJ Z 

X 

X 0 


(XJ 

X 

X 

>r 

X 

4* 

X 

z 

(XJ 

3 vfl 

X 

3 cr cd 

4- 4- 

4* 

4- 

NU 

x cr x x x 4 - 



1 l 



1 3 


3 

< 

1 


3 


1 







1 1 

T V 

Y 

Y 


< 1 





1 


1 



3 



1 

1 













3 

i cm x 

x 

cr 

re qo 

h m 

O 

x 3 00 cr 

X 

X X 

3 

sr 

3 

>4- 


4- 

X 

3 


X 

0 X 

4“ 

X X O' 

3 X 

O 

O' 

vO 

X 0 3 >3- 3 sr 4 - 

3 -T (M 

(XI 

in m co 

x x 0 

(xj <xj x cr 

,— 1 

(XI X 


3 


(XJ 

3 

cr 

(XI 

X 

X 

3 

0 X 

X 

3 X O' 

3 in 

(XJ 

3 

4 - 

3 (xj (\j cr (xj x 

3 X re 

X 

nT >0 m 

>4- X >4- 

x >4” x x 

X 

X X 

a 

<r x 

X 

<r 

X 

X 

X 

(XI X 

X X 

X 

X X 4" 

X 4" 

X 

X 

4- 

XCLXXXX4 

0. 

0. 

0. 

000 

O O 

OOOOOO 

0 0 X 

0 0 

O O 

00000000000 

00000 

oxooooo 










CJ 

















0 

4 00 ^ 

• • • 

4 


X 3 

O' 

fXI 

cr x x cr 

X 

3 

6 


O' 

O' 



cr 

3 <r 

X 

X 

cr 4 - 

4- 

V 

8 

9 

8 

8 

3 

r-i 

X 

X x x (XJ 0 0 


* 


• 

• | • « 

• 

« « 

O 

• 

• 


• 

• 

• 

9 

♦ 

• 

♦ • 

• 

• • • 

• • 

• 

• 

• 

• 0 • • • • • 

' 3 c\J 

(NJ 


x x 0 

rv O O x 

X 3 3 


O' 

X 

O GO 

X 

X 

O' 

X 

CO 

(x> cr 

3 

x 3 in 

4- 3 

(XJ 

(XJ 

■~o 

4 - OD X X X x 

r <n m 

no 


>r -r 

nO (M r- 

- 4 - <xi cu r-* 


X X 


X 

r M 

X 

X 

X 

(XJ 

4- 

X 


1 X 

X 

4- CM (XI 

X X 

X 

X 

(XJ 

(XI X CXI X X x 

-1 ^ 






l 3 1 






| 

1 






1 



1 1 

1 

1 






1 1 





1 

O 





1 







1 1 





C 1 1 










z 

















z 

r> 3 -1 

O' 

cr 

sT r-4 

PXI O' 

>4- 

X x 0 cr 

cr 

x 


X 

X 

— 

3 

O' 

0 

fXJ 

(XJ 

cr 

X O 3 

3 4- M 

<XJ 4- 

4 

4- 

4- 

X (XI X 3 X 3 

ro®s 

fXJ 

in in 

s}- -< 

X 

<xj x x cr 

0 

3 <T 

* 

c* 

X 

•4* 

(X) 

X 

X 

(XJ 

O' 

X 

X X 

cr 

4- X 3 

P^ O 

(4 

(XI 

4) 

4- * 3 0 X (XI X 

j in ro 

no 


so in 

fXI (XJ 


x i r- 

(\| 

<r 3 


X 

X 


(XI 

X 

(XJ 

X 

1 

(XJ 


4" 

X (XJ 4" 

4" fxi 



X 

re (XI (XJ X X 4" 


xx 4 'XCrx 3 fxjxxcr 3 cNj(\j 4 'XX(xi 34 -xxxx<Ni 3 xx 3 crx 3 xxcrr\j(Nj 30 P'-QcooxNr 

larger— 1 — txc'cxjxcxjaocrx 




, _i ^ -r 4 ■ 




>Offioo®xo®>OHioco(fla'«o^Mr\( i ninom^'i'H 3 vM 5 '(M(T'M^h ^ 1 
3 | I 3 k-« (\| ^ I - — • I f— • | I MCOM I I 1 t - 1 CM - ‘ 

I II II I 


I 3 


imO'^'CvD'OtNJ'OX 

* *— • i i x i I rsj 
l i 


z 

o 


'0C0NN07''toain>0^rn'0M^'t7'^O'jfvj'na''4 , Ncr>i-ino^'t®>0(M-i'0^®®f< , SLn>0(«^ 

(\imxx''nfl03crxx(Nj>rx3(.r.-*3crxr'xaoo343'\jx(xjoP , ‘-oo*-(XX('nxiNjX3XOOX 

mr~' 0 ' 0 <jNr '-4 ^ w ^ r- i 4 m xenr— x<xi43xx i invrncos 1 0 'Oin>j(MtNU'ivO‘b^^r N '0 
I III 1 II I I I I I I I 

■-•r“r s 'hO(7''0®l7'f\lC'roj'^''4''l‘0''0 ( Mr^H^O- l < r lf > -lf'^NOCiO'dUr^<-i’J30>-ir'J'4'T'f r '*CPi r '(7' 
r o N j-r-3X03X3nox— <30'3>i'3'D3X4X3XOX(Nja'3XOXr\j404 44xx30xcox 

•Tifi'f^ro»oinr*iir\4-'T^io(ti'}' s j->rioin^c v jr^ir\stLn(M'rrO(Min'£)ir'in'j'r^N(Mf r isf'ii - 4 - x x x x 

0000 0 000000000000000 0 000000000000000 3 00000000 

xoxinoo30ocxoo:coxoxrNjoxr\irNjOOX30Xc\j33xxcooxoOxXL' s X(NjxoG 

3 X 4 ir(xjr- 3 tr 3 croooo'X>(xi'^ 0 '''~ 4 (Nj(N.xxx — iM(MN' 0 'Of'' f *cr' 030 ^^ir' 0 \*Nr 

i *- H '0'0 r '~cr'4-^'X)ocror\j'7Ncn<3'0ooO'-^^nor-O'>0 , - Hr ~^0f^-rsiO<r 

| rxj P — r — | vO in O' N -t sO nj r-< m t — iNvOONO'OvC'I' 4 X O ^ !\l (\ T s - 

3 3 3 3 | ,— < — i i — » • i ii- 1 i it ^ ^ i i i i ^ ^ ^ 7 1 ^ 

II II II 1 

oao3aooorxjX3aooC(\;o(xJOOXOGocoxoaDX03X3cox<Njr\jX(\jxoxxxoxxo 

rN.moC'oofsjp'— O'-iinuO'iruif'^J'O' - 1 ^ t vO l'uu md ^ ci in in ^ o r 'T ° a ' * ui 3 <r 

u' -T u ip >o f 1 ' >11 nuf 1 u u cr o cu >x cu 1 — ' <r <r r~- og • — • uu ■ — i n- i^- 'O ii.’ m f\i (M tr ct i ^ 1 X 1 r^-i_; 

I o , nv!Min^in'nmr'-( j- mHfviSvu^njLn'T --r 43ce>(NjmxX4X3<r(Nj(Nj3>r3XX4ru"' 
III t II II I I I I I 


an z 

l_j LJ 

a. o 

>— <L 
3 > I 

oc a <t 

o cc cc 

CL CL d 


O < < 

Z cC uj 
O Q O 
XHZ 
o UJ — 

cc cc cc 
cl cl a. 


z a . 

<i 

x. < \ 


I U l/V 

■ o £-J 

■ LU <' r 


■ c uj 

> >• 3 
. CL Q 


<3 (_> 

L Z 

a o ; 
■Zkuj 
■ UJ z cc 

I UJ UJ LU 
3 3 3 

□ O Q 


U 

_| < UJ 
mCUI 

3 c <r 
a oc oc 


O U M l 

Z 3 > : 
U 3 -5 i 

• 3 _J O . 

a a * 

Of l/> > 
< UJ LU 

a: a: a: 


Z “vl Z 


a ; 


<r O 


3 f 
O 


o c 

LU *-H 

oc qc 


OnczocuarO Z z < _ 

3 cc 3 o cc o o <t l: a uj «g 3 x ^ 
<r > uu <X3CCX00<0 

03 a<u.<M^<Q<'nni 

<1 < 3 3 L— L- LU iC Cl CC 
O > X >^ZZZZ3L0CLOZZOCJ 
►-iOOCJ<r<L<<<j<<<<i<<i< 
cc cc lai/ii/ii/ii/ii/ii/ii/M/ii/M/iwuo 


x cc 


« x 3 

J X o 

: o o 
> 00 x 


<3 3 3 a. 
o 3 t — | — 1 
z u «i x 


26 


SE4GR0VE -37 . Cti 


Table 2 

Tracings from Surface to Surface 
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Table 2 

Tracings from Surface to Surface 
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Table 3 

Tracings 3Q0 Km Above Stations 
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Table 3 

Tracings 300 Km Above Stations 
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Table 3 

Tracings 300 Km Above Stations 
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Table 4 

Tracings 400 Km Above Stations 
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Tracings 400 Km Above Stations 
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Table 4 

Tracings 400 Km Above Stations 
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Table 5 

Conjugate Pairs (with 1.1 < L < 12.) 
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Table 5 (Continued) 
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onolulu - Noumea 1.15 2.3 0.60 



Table 6 

Stations with Conjugate Areas in the Antarctica (L < 12.) 


Station 

L 

Conjugate Area 

Kap Tobin 

10.94 

Enderby Land 

Barter Island 

9.05 

Ocean Station (South Pacific Ocean) 

Fort Chimo 

8.40 

South Polar Plateau 

Cape Jones 

6.58 

Near Byrd Station 

Knob Lake 

6.35 

Edith Ronne Land 

Narssarssuaq 

7.50 

Queen Maud Land 

Julianehaab 

7.31 

Queen Maud Land 

Goose Bay 

5.33 

Edith Ronne Land 

Mont Joli 

4.26 

Ellsworth Highland 

Val D'Or 

4.22 

Mary Byrd Land 

Caribou 

3.86 

Near Eights Station 

Deep River 

3.78 

Mary Byrd Land 

Fort Monmouth 

2.82 

Peter I Island 
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Table 7 

Conjugate Areas for Stations with 1.1 < L < 12. 
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Corinaldo 1.65 Boshof (Orange Free State) RSA 

Cook 1.74 Birakan (near Birobidzhan) USSR 

Curacao 1.27 San Carlos de Bariloche Argentina 
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1.14 Auvergne (Northern Territory) Australia 


Table 7 (Continued) 
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Nitsanim 1.18 Near Vila Cabral Mozambique 

Nouvelle Amsterdam 2.32 Demyanovka (Kazakhstan) USSR 

Onagawa 1.31 Butru (Queensland) Australia 

Oohira 1.23 Near Camooweal (Queensland) Australia 

Paramaribo 1.18 Near Florianopolis Brazil 
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Table 7 (Continued) 
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LINES OF CONSTANT L (earth radii) 
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LINES OF CONSTANT L (earth radii) 
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LINES OF CONSTANT L (earth radii) 
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LINES OF CONSTANT L (earth radii) 



ALTITUDE 2500 KM 
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LINES OF CONSTANT L (earth radii) 
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LINES OF CONSTANT L (earth radii) 
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LINES OF CONSTANT B (gauss) 
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LINES OF CONSTANT B (gauss) 
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t 


ILII 


IB&11I 


M 




m si iwim 

mmiig 

mmmamvmmm 
wmmniMmwimm 

■UBS! 


imwimmwA 


I 'BLfl 


ISSIHRII 

Bill'll 


iiiiihb 

lllllklB 


mmm 


■RMlBBBniWIIIB^^IIiniKa^HHBiaBHl 
HBUHKHIHniRi«B«MBBBnSIHHiBi 

MiMEMBi MBiBigSMB 

KnB^ESnViHBjBHHKflHSuBHfiSilESl 

^mwanmummmmumumnmumuusmusmumsmuKmmsim^a 

|BHnHHn| 
IBHK«EHSSniHHflKHnHttSillSSkBE3 

ligi^iMBigMB MMgaBrji BwnBKsagagBSgi i mi ■■■ i ■iMg 

sail il%r4 


b gw jMBiinn| 

HHHIsaBiBBHHHHU 


iii~~ Bsa Bffli 

I jHiw" 

pj« 
leu? 


ikbb 


■BUBBIHBMIiiiaBnHmilTOBIHUMi 
■BlKRSBnimHftlinnHIliMBHninilBIBBBlHBi 

mmmw^mamamuamkmmaimmmmMmmwnnmummmiamuml 
1 1 — m mm bbw ■■ ihhibbbbbi bbiimwbbmibbi I 

BBBBBBI 8 imBI BmMMHl 
MjgBBgBMBBBIMMM 

|lBMMB«B?&£flBBBMBL11*iklimiBBr ... 


vmamk 


m 


^gg 

H 

■HBfflBMKBHBIIflBBHlIM 

grjj 

IB 

n 


MB 

H 

H 

BKI 

IBB 

§H 


mmmrsmummuwmmmuuunnmu 

KB 

K 

■ 

B 


HUS 

IB 

H 


^BH 


mismM’mmmummmmmmuuummwm 

mb 

rj 

IK 

HI 

u 

H 

WM 

H 


H 

SI 


■iMBWiavniiiimimrii 

ins 

ini 

B 

BH 

B 


■i 

HE 




ALTITUDE 1400 KM 






LINES OF CONSTANT B (gauss) 



ALTITUDE 1500 KM 




LINES OF CONSTANT B (gauss) 




ALTITUDE 1600 KM 






LINES OF CONSTANT B (gauss) 



ALTITUDE 1700 KM 


LINES OF CONSTANT B (gauss) 
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ALTITUDE 1800 KM 






LINES OF CONSTANT B (gauss) 



ALTITUDE 1900 KM 




LINES OF CONSTANT B (gauss) 



ALTITUDE 2000 KM 








LINES OF CONSTANT B (gauss) 



ALTITUDE 2100 KM 




LINES OF CONSTANT B (gauss) 
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ALTITUDE 2200 KM 






LINES OF CONSTANT B (gauss) 



ALTITUDE 2300 KM 





LINES OF CONSTANT B (gauss) 
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ALTITUDE 2400 KM 











LINES OF CONSTANT B (gauss) 
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ALTITUDE 2500 KM 


LINES OF CONSTANT B (gauss) 



ALTITUDE 2600 KM 




LINES OF CONSTANT B (gauss) 
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ALTITUDE 2700 KM 






LINES OF CONSTANT B (gauss) 



ALTITUDE 2800 KM 





LINES OF CONSTANT B (gauss) 
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ALTITUDE 2900 KM 






LINES OF CONSTANT B (gauss) 



ALTITUDE 3000 



